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ABSTRACT 


Winter  1976  Data  Systems  Test  wind  data  was  used  to  produce  a 
number  of  global  kinematic  flel^^  The  data  was  prepared  by  the 
Goddard  InstUiite  ^Por-Space  Studies  Atmospheric  General  Circulation 
Model.  Interpolation  scheme  was  developed  to  move  the  wind  data 
from  sigma  levels  to  constant  pressure  levels.  Time-averaged  fields 
of  the  zonal  wind  component,  meridional  wind  component,  divergence, 
vorticity,  steady  energy  and  transient  energy  were  computed  for 
200  mb,  500  mb  and  850  mb.  In  addition,  the  standing  kinetic  energy 
and  transient  kinetic  energy  were  decomposed  by  wavenumbers  Into 
energy  spectra  for  200  mb  and  500  mb  for  different  latitudes. 

Results  of  this  research  were  compared  with  other  similar  studies 
and  showed  favorable  correlation.  Resolution  of  synoptic  scale 
features  was  good.  Organized  areas  were  found  In  the  time-averaged 
divergence  fields  Indicating  that  divergence  features  can  be  re- 
solved In  this  data.  This  study  shows  that  the  Goddard  Model  Is 
well  suited  for  globakstudles  of  the  general  atmospheric  circulation. 
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CHAPTER  I 
INTRODUCTION 

The  general  circulation  of  the  atmosphere  is  still  not  well  under- 
stood. Lorenz  (1967)  mentions  that  even  though  there  have  been  over  a 
thousand  excellent  studies  on  the  subject,  a full  explanation  of  the 
distribution  of  the  atmospheric  wind  patterns  is  lacking.  The  atmos- 
phere is  so  vast  and  motions  in  it  are  so  complex  that  it  was  impossible 
to  collect  global  circulation  data  until  this  century.  It  is  obvious 
from  our  current  knowledge  of  the  atmosphere  that  to  successfully  un- 
ravel this  mystery  we  need  more  clues. 

The  clues  in  this  case  are  more  complete  data  sets.  Data  is  need- 
ed for  the  entire  globe  at  several  levels  in  the  atmosphere.  The  only 
way  that  the  required  volumes  of  data  can  be  collected  is  through  inter- 
national cooperation  of  countries.  Although  much  data  is  collected  and 
exchanged  daily  by  many  countries,  it  is  not  sufficient  for  the  task  of 
fully  explaining  the  general  circulation. 

Meteorologists  recognized  several  years  ago  that  it  was  time  to 
formulate  a plan  for  gathering  more  detailed  global  data.  As  a result, 
the  Global  Atmospheric  Research  Project  (6ARP)  was  formulated.  GARP 
will  be  the  most  intensive  effort  yet  to  collect  and  analyze  atmospheric 
data  for  the  entire  globe.  It  will  involve  many  countries,  thousands  of 
scientists  and  much  sophisticated  equipment  when  it  is  implemented  at 
the  end  of  this  decade. 


To  be  useful,  the  massive  amount  of  data  which  will  be  generated 
by  GARP  will  have  to  be  placed  into  a format  which  will  allow  computer 
analysis.  The  Goddard  Institute  for  Space  Studies  (GISS)  Nine-Level 
Atmospheric  General  Circulation  Model  has  been  selected  for  this  data 
preparation.  This  study  uses  a set  of  data  generated  by  the  GISS 
model . 

The  objective  of  this  study  was  to  compare  the  data  prepared  by 
the  GISS  model  to  data  from  other  sources.  This  was  done  partly  to 
assess  the  quality  of  the  data  produced  by  the  model  and  partly  to  get 
a feel  for  what  type  features  can  be  resolved  by  this  data.  The  anal- 
ysis of  the  data  was  divided  into  two  basic  areas.  First,  a number  of 
time-averaged  fields  were  generated.  These  included  time-averaged  wind 
components,  divergences  and  vorticities.  The  second  area  of  study  in- 
volved a wave  number  decomposition  of  the  kinetic  energy  for  the  stand- 
ing and  transient  eddies.  The  period  covered  by  the  data  (one  month) 
was  too  short  to  consider  any  attempt  to  add  to  the  theory  of  the 
general  circulation.  However,  this  study  gives  an  indication  of  the 
type  and  amount  of  data  which  will  soon  be  available. 


Global,  synoptic  scale  wind  data,  which  covers  35  days  at  12  hour 
intervals,  was  used  for  this  analysis.  The  data  covers  the  period  from 
OOZ/29  JAN  76  until  12Z/3  MAR  76.  It  was  prepared  by  the  GISS  atmos- 
pheric circulation  model  and  is  designated  Data  System  Test  (DST)  data. 
The  model  uses  five  primary  quantities  to  represent  the  behaviour  of 
the  atmosphere:  u—the  zonal  wind  speed  (west  to  east),  v— the  merid- 
ional wind  speed  (south  to  north),  T“the  temperature,  q—the  specific 
humidity  of  water  vapor,  and  w— the  surface  pressure  minus  the  pressure 
at  the  top  of  the  model.  This  study  used  only  the  u,  v,  and  ir  data. 
Three  forcing  terms  are  included  in  the  differential  equations  used  by 
the  model  to  update  the  primary  quantities.  These  terms  are  designed 
to  represent  the  physical  processes  of  earth-atmosphere  and  sun-atmos- 
phere interaction  and  include  a horizontal  frictional  force  term,  a 
diabatic  heating  rate  term,  and  a rate  of  moisture  addition  term.  For 
detailed  information  on  the  GISS  model  consult  Somerville  et  al.  (1974) 
The  model  uses  a split  global  grid  in  spherical  coordinates.  The 
grid  is  diagramed  in  Figure  1.  The  primary  grid  points,  denoted  by 
o's,  contain  T,  q,  and  ir  data  while  u and  v data  is  carried  at  the 
secondary  grid  points,  the  +'s.  The  poles  are  primary  grid  points, 
while  the  equator  is  at  secondary  grid  points.  Grid  points  are  spaced 
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both  grids.  This  makes  a total  of  3312  points  over  the  globe  per  grid. 
The  primary  grid  has  46  points  from  pole  to  pole  and  72  points  around 
a latitude  circle. 

The  secondary  grid  contains  the  same  number  of  points,  but  the  row 
which  falls  below  the  South  Pole  is  meaningless  so  it  only  has  45 
meaningful  latitude  circles.  Figures  2 and  3 show  the  primary  grid 
for  the  Northern  Hemisphere  and  Southern  Hemisphere  respectively  with 
the  topography  used  by  the  model  given  in  meters. 

A sigma  vertical  coordinate  system  is  used.  Sigma  ( a ) is  de- 
fined by 


where  P is  the  pressure  at  the  level  of  measurement, 

Pg  is  the  surface  pressure  at  the  grid  point, 

Pj  is  the  pressure  at  the  top  of  the  atmosphere. 

Note:  P-j-  is  10  mb  in  this  model.  There  are  nine  levels  above  the  sur- 
face with  sigma  at  each  level  (o^)  given  by  Equation  (2.2). 

^ t = 1,  2---9.  (2.2) 

Level  one  is  the  highest  level  and  nine  is  the  lowest.  Note  that  a is 
zero  at  the  top  of  the  atmosphere  and  one  at  the  surface. 

Sigma  Coordinates  Over  Mountains 
Divergence  fields  at  constant  pressure  were  computed  at  several 
sigma  levels  using  the  following  equation  for  divergence: 
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Figure  2 GISS  Model  Northern  Hemisphere  Topography  Terrain  heights 
are  given  In  meters. 
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a r 1 + iv  f 1 3in  i 

IT  '■3a  '•a  cos  BX^  3a  '•a 


(2.3) 


The  finite  difference  form  of  Equation  (2.3)  which  was  used  is 


(V..  V). , . i r— !—  ~ ■ ''i.j-i 

P M,j  a kos  (jij  AX  A(J) 


■ ''i.j  *i~  f 


afl  - uf-' 


j.  r“<.J  ' f 1 "W-J  ■ *1-1  .J' 

(■a  cos  iij  A>  J 


Aa  '•a  cos 


vf!  - v‘-' 


t M.j  ~ M.J  fl  Cl.J^l  ■ 

Aa  '•a  A4>  ^-1  • 


(2.4) 


After  examining  these  fields,  it  was  suspected  that  the  values  over  the 
mountains  were  spurious.  Thus,  the  following  test  was  conducted.  A 
vertical  wind  profile,  constant  at  a given  pressure,  was  selected  and 
entered  at  all  grid  points.  Figure  4 illustrates  the  wind  profile, 

The  divergence  was  then  computed  for  several  sigma  levels  for  the  ir 
field  of  002/8  FEB  76  using  Equation  (2.4).  Values  over  the  oceans  were 
small  (order  10"®  sec"^  or  smaller)  as  they  should  have  been.  However, 
over  mountains  such  as  the  Andes,  values  on  the  order  of  10"®  sec"^ 
were  obtained.  These  values  occurred  because  of  the  steep  slopes  of 
the  sigma  surfaces  as  they  cross  over  the  mountains.  It  is  not  possible 
to  correctly  approximate  the  derivatives  involved  by  finite  differences 
for  the  coarse  grid  used  in  this  version  of  the  GISS  model.  At  this 
point  the  decision  was  made  to  Interpolate  the  wind  data  to  constant 
pressure  levels  to  partially  circumvent  this  problem. 
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Interpolation  Scheme 

Although  the  advantages  of  sigma  coordinates  are  well  known, 
(Phillips, 1957) , having  the  data  on  constant  sigma  surfaces  proved  to 
be  inconvenient  for  this  study.  As  seen  in  the  previous  section,  it 
was  not  possible  to  produce  meaningful  divergence  fields  on  constant 
sigma  surfaces  in  the  vicinity  of  mountains.  Also,  if  the  data  was 
analyzed  at  constant  sigma  levels,  it  could  not  be  readily  compared 
with  other  data  since  most  analyses  of  the  type  performed  here  are  done 
at  constant  pressure  levels.  For  these  reasons  the  decision  was  made 
to  interpolate  the  data  from  the  sigma  levels  to  constant  pressure 
levels.  Three  levels  were  chosen  for  study,  200  mb,  500  mb,  and  850  mb. 

A cubic  spline  was  used  for  the  interpolation  from  sigma  levels  to 
constant  pressure  levels.  The  spline  was  essentially  the  one  contained 
in  Conte  and  de  Boor  (1965).  The  procedure  involves  constructing  a 
piecewise-cubic  polynomial  which  interpolates  the  function  at  each  of 
the  known  points.  This  is  done  by  constructing  a certain  cubic  poly- 
nomial of  the  form 

P,(x)  = ’i.i  <2.5) 

for  each  sub-interval  so  that  the  resulting  piece-wise  cubic  inter- 
polant  is  twice  continuously  differentiable, 

pj'^(x^)  = Pj"(x^)  , i=2,  3,-*.N-l 

where  N is  the  number  of  known  points.  A piece-wise  cubic  interpolant 
of  this  form  is  known  as  a cubic  spline. 

The  subroutine  developed  for  use  in  this  interpolation  required 


n 


knowledge  of  the  function  values  at  each  known  point  plus  the  slope  of 
the  function  at  the  two  end  points.  The  end  point  slopes  were  found  by 
finding  the  derivative  of  a quadratic  polynomial  passing  through  the 
end  point  and  the  two  adjacent  points.  The  subroutine  then  used  Gauss 
elimination  to  solve  for  the  coeffeclents  a^  r 1’  1 ^4  1 

Equation  (2.5).  A minimum  of  four  points  are  needed  for  the  cubic 
spline  of  this  subroutine. 

To  Interpolate  the  winds  to  constant  pressure  levels  It  was  neces- 
sary to  have  ir  values  at  the  secondary  grid  points.  These  values  were 
derived  by  summing  the  ir  values  at  the  four  primary  points  surrounding 
a secondary  point  and  dividing  by  four.  Though  simple,  this  procedure 
produced  acceptable  results. 

Once  the  ir  value  at  a secondary  grid  point  was  known,  the  pressure 
at  each  sigma  level  was  computed.  The  nine  pressure  values  were  then 
passed  to  the  spline  as  the  values  and  the  wind  components  at  each 
level  were  passed  as  the  fCxi)  values.  The  cubic  spline  then  computed 
f(x)  values  at  200  mb,  500  mb,  and  850  mb. 

For  levels  which  were  below  ground,  a value  of  999999  was  assigned 
as  a flag.  The  effect  of  this  flagging  appears  In  the  850  mb  charts 
presented  In  the  following  chapter. 

The  cubic  spline  was  tested  on  several  known  functions.  Table  1 
shows  the  functions,  the  exact  value  of  the  function,  and  the  Interpo- 
lated value  of  the  function,  and  the  error.  From  the  results  one  can 
see  that  the  cubic  spline  Is  sufficiently  accurate  for  any  reasonably 
smooth  wind  profile.  The  results  presented  In  the  next  chapters  sug- 
gest that  this  Interpolation  scheme  was  very  successful. 
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CHAPTER  III 
TIME-AVERAGED  FIELDS 

Time -averaged  fields  of  several  quantities  were  constructed  for 
200  mb,  500  mb,  and  850  mb  levels.  These  fields  are  presented  in  this 
chapter  with  a discussion  of  how  they  compare  with  similar  fields  com- 
puted from  other  data.  It  should  be  noted  that  for  both  hemispheres 
positive  values  of  u denote  flows  from  west  to  east  and  positive  v 
values  indicate  flows  from  south  to  north. 


Wind  Fields 

For  each  grid  point  the  time-averaged  zonal  and  meridional  wind 
components  were  defined  in  the  usual  manner  with  Equations  (3.1)  and 
C3.2). 


T 

/ udt 
0 


(3.1) 


(3.2) 


u = zonal  wind 

u = time-averaged  zonal  wind 

V = meridional  wind 

V = time-averaged  meridional  wind 
T = time  interval 


The  quantities  u and  v were  evaluated  by  Equation  (3.3); 

1 N 

X = N (3.3) 

with  X = u or  V.  For  this  sample  N = 70.  Figures  5,  6,  and  7 
show  the  IT  fields  for  200  mb,  500  mb,  and  850  mb  respectively.  Figures 
8,  9,  and  10  display  the  corresponding  v fields.  The  asterisks 
on  the  850  mb  charts  represent  grid  points  where  the  850  mb  surface  Is 
below  the  surface  of  the  earth  In  the  model. 

At  each  grid  point  the  standard  deviation  (SO)  for  the  above 
fields  was  also  computed  as  follows: 


where  N = 70  and  x^  = u^-  or  v^. 

Figures  11  and  12  present  the  standard  deviations  of  the  IT  and  v 
fields  respectively  at  500  mb.  The  standard  deviations  at  the  other 
levels  showed  similar  patterns  so  they  are  not  Included. 

A comparison  of  the  time-averaged  fields  with  fields  from  other 
data  sets,  (Van  Loon,  1971),  (Sadler,  1975),  showed  good  agreement  and 
Indicated  that  the  time-averaged  wind  patterns  for  February  1976  con- 
formed closely  to  expected  values.  Some  points  of  Interest  on  the  iT 
charts,  include  the  band  of  pronounced  westerly  flow  at  all  levels  In 
the  Southern  Hemisphere  centered  around  48°  S latitude,  the  strong 
westerly  flow  at  200  mb  In  the  Pacific  centered  about  4°  S,  and  the 
easterly  flow  at  200  mb  in  the  Southern  Hemisphere  tropics  which  seems 
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Figure  6 Time-average  of  u wind  component  at  500  mb  (m/sec) 
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Figure  7 Time-average  of  u wind  component  at  850  mb  (m/sec); 

asterisks  indicate  areas  where  850  mb  surface  is  below 
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Figure  9 Time-average  of  v wind  component  at  500  mb  (m/sec). 
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Figure  10  Time-average  of  v wind  component  at  850  mb  (m/sec); 

asterisks  indicate  areas  where  850  mb  surface  is  below 
ground. 
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Figure  11  Standard  deviation  of  time-averaged  u wind  component  at 
500  mb  (m/sec). 
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Figure  12  Standard  deviation  of  time -averaged  v wind  component  at 
500  mb  (m/sec). 


to  correspond  with  the  continents.  Also,  the  u fields  show  that  the 
expected  maximums  of  westerly  flow  in  the  Northern  Hemisphere  are 
visible  at  all  levels.  They  occur  along  the  eastern  United  States, 
over  northern  Africa,  and  off  the  coast  of  Asia. 

Extensive  easterly  flow  is  evident  at  850  mb  throughout  much  of  the 
tropics.  This  850  mb  easterly  flow  shows  a distinct  maximum  near  the 
equator  under  the  maximum  of  the  200  mb  westerly  jet.  This  may  repre- 
sent the  Walker  circulation  described  by  Krishnamurti  et  al.  (1973). 

At  all  levels  standard  deviations  of  IT  are  generally  smaller  than 
the  corresponding  value  of  u.  For  the  200  mb  and  500  mb  levels  the  SO 
values  are  25  - 50  % of  IT.  The  standard  deviations  of  V are  generally 
larger  than  v at  all  levels.  Many  SO  values  are  twice  the  value  of  7 
and  sometimes  the  SO  is  an  order  of  magnitude  larger  than  7.  The  fact 
that  V shows  more  variability  than  IT  is  consistant  with  a pattern  of 
waves  traveling  around  the  latitude  circles. 

The  V fields  show  a high  degree  of  vertical  continuity  when  the 
three  levels  are  compared  and  indicate  good  resolution  of  the  vertical 
characteristics  of  the  waves.  A comparison  of  the  long  wave  patterns 
with  the  continents  supports  the  hypothesis  that  large  land  masses  play 
a primary  role  in  the  production  of  the  long  waves. 

It  is  interesting  to  note  that  there  are  many  climatological  charts 
of  zonal  flow  available  for  many  levels,  but  global  charts  of  the  time- 
averaged  meridional  flow  are  not  readily  available.  More  charts  of  this 
type  will  probably  appear  because  of  the  relative  ease  of  producing  them 
from  a global  model  such  as  this. 

Longitude-and-time  averaged  values  of  u and  v,  which  are  defined  by 
Equation  (3.5)  were  computed  using  Equation  (3.6). 
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f2n  fT 

= 27r  /o  /o  ’“‘TdX 


(3.5) 


1 


M N 


NM  <*t>J 


(3.6) 


where  X = longitude, 

(xi)j  = (ui)j  or  (vi)j. 

M = 70  (number  of  time  periods) 
and  N = 72  (number  of  longitudinal  grid  points). 

Figure  13  is  a vertical  cross  section  for  Uq.  Values  were  comput- 
ed at  200  mb,  500  mb,  and  850  mb.  The  x's  along  these  three  axes 
represent  the  points  where  the  isotach  values  occurred.  The  M's  along 
the  axes  indicate  points  where  relative  positive  maximums  occurred,  and 
the  N's  indicate  where  the  largest  absolute  values  of  negative  Uq's 
occurred.  The  results  for  Vq  at  the  three  levels  are  tabulated  in 
Table  2. 

Figure  13  agrees  very  closely  with  results  obtained  by  Heastie 
and  Stephenson  (1960)  for  January  1958.  There  is  a slight  indication 
of  both  the  polar  jet  and  the  subtropical  jet  at  500  mb  in  the  Northern 
Hemisphere.  In  the  Southern  Hemisphere  the  only  jet  indicated  is 
located  about  48“S  latitude. 

The  Vq  values  are  generally  smaller  than  those  of  Sadler,  (1975), 
and  Oort  and  Rasmusson  (1971),  for  the  tropical  latitudes,  but  are 
comparable  to  those  which  may  be  inferred  from  Fig.  1 of  Krishnamurti 
et  al,  (1973).  There  is  good  indication  of  the  northern  Hadley  cell 
in  the  Vq  field  and  also  of  the  position  of  the  ITCZ  at  about  20“S 


Ffgure  13  Vertfcal  cross  sectfon  of  time-and-longitude  averaged  u wtnd  component  (m/sec) 
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at  850  mb.  It  also  appears  that  an  equator-ward  slope  with  height  is 
detectable  for  the  ITCZ. 

Time-Averaged  Kinetic  Energy  Fields 
The  standing  kinetic  energy  per  unit  mass  (SKE)  at  each  point  is 
defined  as: 

—2  —2 

SKE  = . (3.8) 

SKE  fields  for  the  three  levels  were  computed  from  the  time-averaged 
wind  fields.  Figures  14  and  15  are  the  200  mb  SKE  fields  for  the 
Northern  and  Southern  Hemispheres  respectively.  Figures  16,  17, 

18  and  1 9 are  the  corresponding  fields  for  500  mb  and  850  mb. 

These  charts  show  the  average  jet  stream  positions.  Three  dis- 
tinct jets  are  detectable  in  the  Northern  Hemisphere  in  mid-latitudes 
at  200  mb  and  500  mb.  One  along  the  East  Coast  of  the  United  States, 
one  over  northern  African  and  the  Middle  East,  and  the  strongest  one 
off  the  coast  of  Asia.  At  850  mb  there  are  only  two  wind  maxima,  the 
one  off  the  East  Coast  of  the  United  States,  which  is  the  strongest, 
and  the  one  off  the  coast  of  Asia.  The  500  mb  chart  shows  some  indica- 
tion of  the  flow  splitting  and  going  around  the  Tibetian  high  country. 

The  locations  of  the  jet  cores  are  in  good  agreement  with  climatologi- 
cal charts  for  January  and  February,  (Namias,  1963),  (Sadler,  1975). 

In  the  Southern  Hemisphere  only  one  mid-latitude  jet  is  detectable 
but  it  appears  at  all  three  levels.  This  jet  is  centered  along  the 
48“  S latitude  line  and  runs  west  to  east.  This  is  in  constrast  to  the 
Northern  Hemisphere  jets  which  are  oriented  from  southwest  to  northeast. 

The  position  of  this  jet  core,  centered  about  48“  S latitude  and  between  5 
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Figure  14  Standing  kinetic  energy  per  unit  mass  for  N.H,  at  200  mb 
(mVsec^);  dashed  isolines  are  for  values  of  500,  1000, 
1500  and  2000;  dotted  isoline  denotes  a value  of  750. 
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Figure  15  Standing  kinetic  energy  per  unit  mass  for  S.H.  at  200  mb 
(m^/sec^);  dashed  isolines  are  for  a value  of  500,  dotted 
isolines  denote  750  values. 
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Figure  16  Standing  kinetic  energy  per  unit  mass  for  N.H.  at  500  mb 
(m^/sec^);  dashed  isolines  are  for  values  of  200,  300, 
400  and  500. 
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Figure  17  Standing  kinetic  energy  per  unit  mass  for  S.H.  at  500  mb 
(m2/sec2);  dashed  isolines  are  for  values  of  200,  300 
and  400. 
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Figure  19  Standing  kinetic  energy  per  unit  mass  for  S.H.  at  850  mb 
(m^/sec^);  dashed  isolines  are  for  values  of  50,  75,  and 
100;  areas  where  850  mb  surface  is  below  ground  are  blank. 


37“  E and  78“  E longitude  is  also  in  agreement  with  other  similar 
analyses,  (van  Loon,  1971). 

In  the  tropical  Pacific  area,  the  westerly  jet  at  200  mb,  which 
crosses  the  equator  is  of  interest.  The  core  is  at  4“  S latitude  and 
the  wind  speeds  are  considerably  stronger  than  the  February  average 
given  by  Sadler,  (1975).  He  shows  about  45  knots  as  the  maximum  where 
this  data  has  an  area  with  wind  speeds  of  near  60  knots  with  the  maxi- 
mum about  65  knots.  An  examination  of  the  IT,  v and  steady  energy 
charts  for  200  mb  reveals  a large  Pacific  trough  associated  with  this 
jet. 

Equation  (3.8)  defines  the  time-averaged  transient  kinetic  energy 
(TKE)  per  unit  mass  as 


TKE  = i / t-—  dt  (3.8) 

' 0 ^ 

with  u*  = u - u and  v*  = v - v. 

The  TKE  was  computed  directly  from  the  standard  deviation  fields 
of  the  time-averaged  wind  components  as  follows: 

o , N , 

^ = F I (x-  - x)^  (3.9) 

N i*i  1 

Xj  = Uj  or  v^  and  N = 70. 

TKE  = SD.U_^  S^.  (3  IQj 

where  SOU  and  SDV  are  the  standard  deviations  of  u"  and  v respectively. 
These  fields  are  depicted  in  Figures  20,  21,  22,  23,  24  and  25, 
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Figure  21 


Time-averaged  transient  kinetic  energy  per  unit  mass  for 
S.H.  at  200  mb  (m^/sec^);  isolines  are  for  SKE  and  are 
from  Figure  15. 
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Figure  22  Time-averaged  transient  kinetic  energy  per  unit  mass  for 
N.H.  at  500  mb  (m^/sec^);  isolines  are  for  SKE  and  are 
from  Figure  16. 
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Figure  24  Time-averaged  transient  kinetic  energy  per  unit  mass  for 
N.H.  at  850  mb  (mVsec^);  isolines  are  for  SKE  and  are 
from  Figure  18;  areas  where  850  mb  surface  is  below  ground 
are  blank. 
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Figure  25  Time-averaged  transient  kinetic  energy  per  unit  mass  for 
S.H.  at  850  mb  (m^/sec^);  isolines  are  for  SKE  and  are 
from  Figure  19;  areas  where  850  mb  surface  is  below  ground 
are  blank. 
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The  isotachs  from  the  corresponding  standing  kinetic  energy  charts 

are  Included  on  the  TKE  fields,  This  makes  it  possible  to  compare  :j 

, j 

the  transient  energy  with  the  average  position  of  the  wind  maxima. 

By  using  the  u charts  it  is  possible  to  compare  the  transient  IJ 

energy  on  the  upwind  and  downwind  sides  of  mountains  in  several  loca-  |] 

’ I 

tions.  By  following  along  the  44°  N latitude  circle  it  is  possible  to 
detect  a drop  in  the  transient  energy  to  the  east  of  the  Rockies  and 
the  Himalayas.  This  suggests  that  the  mountains  are  responsible  for 
deflecting  the  transient  eddies. 


Time-Averaged  Divergence  and 
Relative  Vorticity  Fields 

Divergence  and  relative  vorticity  fields  were  produced  from  the 
time-averaged  wind  fields.  In  spherical  coordinates,  the  horizontal 
divergence  and  vorticity  along  a constant  pressure  surface  are  given  by 


DIV 


1 M + 1 iv  _ V tan  <ti 
a cos  (ji  3X  a d(|i  a 


(3.11) 


and 


wno  - 1 IV  _i  3u  ^.u  tan  ^ 

a cos  <ti  ax  a 8(|)  a 


(3.12) 


where  a is  the  radius  of  the  earth,  X is  longitude  and  (|>  is  latitude. 

To  compute  the  divergence  and  the  vorticity.  Equations  (3.11)  and 
(3.12)  were  put  into  the  following  finite-difference  form; 

a '■  4X  cos  (|>j  A(j.  ''ij 


(3.13) 


42 


VOR 


t.J 


a I ax  cos  ♦j  a^  ^ *jj  • 


(3.14) 


Centered  space  differencing  was  used  to  approximate  the  derivatives. 

To  aid  In  the  Identification  of  clusters,  these  fields  were  out- 
put In  two  parts.  One  part  has  the  negative  values  blanked  out  and 
the  other  omits  the  positive  values.  The  positive,  time-averaged  di- 
vergence fields  are  contained  In  Figures  26  , 27  and  28.  Figures 
29  , 30  and  31  show  the  corresponding  negative  divergences.  Posi- 
tive time-averaged  vorticitles  for  500  mb  are  In  Figure  32,  while 
Figure  33  has  the  negative  values.  The  vorticity  charts  at  200  mb 
and  850  mb  showed  similar  patterns  so  they  were  not  Included.  It 
should  be  noted  that  positive  relative  vorticitles  correspond  to 
cyclonic  vorticity  In  the  Northern  Hemisphere  and  anticyclonic  vortici- 
tles In  the  Southern  Hemisphere.  Standard  deviations  for  the  above 
charts  were  calculated  using  Equation  (3.4)  with  replaced  by  the 
divergence  or  the  relative  vorticity.  Figures  34  and  35  are  the 
standard  deviations  of  the  time-averaged  divergences  and  vorticity 
fields  at  500  mb. 

A comparison  of  the  time-averaged  divergence  fields  with  their 
standard  deviations  shows  that  the  standard  deviations  are  generally 
larger  than  the  time-averaged  divergences.  At  some  points  the  standard 
deviations  are  an  order  of  magnitude  larger.  This  suggests  that  the 
time-averages  are  small  numbers  resulting  from  oscillations  between 
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Figure  27  Positive  time-averaged  divergence  at  500  mb  (x10~®sec"M. 
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Figure  34  Standard  deviation  of  time-averaged  divergence  at  500  mb 
(xlO“°sec“M. 
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Figure  35  Standard  deviation  of  time-averaged  relative  vorticity  at 
500  mb  (xlO"°sec"' ) . 
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relatively  large  positive  and  negative  divergence  values. 

Although  the  time-averaged  divergence  fields  do  not  show  patterns 
which  are  as  organized  as  the  patterns  on  the  time-averaged  vorticity 

f 

charts,  there  are  several  clusters  of  interest.  These  clusters  are  best 
seen  at  200  mb  where  organized  areas  of  positive  devergence  are  present 
in  the  tropics  over  South  America,  Africa  and  in  the  vicinity  of  Aus- 
tralia. These  areas  are  similar,  except  for  an  eastward  shift  of 
about  30"  longitude,  to  those  which  may  be  inferred  from  Fig.  1 of 
Krishnamurti  et  al,  (1973).  A smaller  cluster  is  found  over  and  south 
of  India.  Also,  there  is  some  evidence  of  organized  patterns  associated 
with  the  Northern  Hemisphere  jet  streams  in  the  Pacific  and  Atlantic 
Oceans. 

In  the  Southern  Hemisphere  tropics,  it  is  interesting  to  note  that 
the  largest  standard  deviations  of  the  time-averaged  divergences  gener- 
ally occur  in  conjunction  with  the  clusters  mentioned  above.  GOES 
movie  loops  show  that  these  regions  in  South  America  and  Africa  often 
experience  extensive  diurnal  convective  development.  Thus,  one  may 
speculate  that  the  large  standard  deviations  in  these  areas  reflect 
diurnal  oscillations  associated  with  the  growth  and  decay  of  large  scale 
convection. 

Distinct,  organized  patterns  appear  in  the  time-averaged  relative 
vorticity  charts  in  both  hemispheres.  These  patterns  correlate  well 
with  the  jet  stream  positions  which  produces  confidence  in  the  credi- 
bility of  the  time-averaged  vorticity  fields.  The  standard  deviations 
of  these  fields  are  of  the  same  order  of  magnitude  as  the  fields  them- 
selves. These  patterns  imply  that  the  vorticity  fields  do  not  vary  with 
time  as  much  as  the  divergence  fields. 


CHAPTER  IV 


KINETIC  ENERGY  OF  THE  STANDING  WAVES 


The  distribution  of  the  kinetic  energy  of  the  standing  eddies 
and  the  transient  eddies,  as  a function  of  wave  number  was  computed. 
This  was  done  for  wave  numbers  0 to  36  at  both  200  mb  and  500  mb  by 
using  a Fourier  analysis  to  decompose  the  time-and-longitudinally 
averaged  kinetic  energy  into  an  energy  spectrum.  Contributions  to 
the  kinetic  energy  from  the  zonal  and  meridional  wind  components  were 
calculated  separately. 

The  energy  spectrum  for  the  transient  energy  is  presented  in 
Chapter  V,  while  the  spectrum  for  the  standing  eddies  is  covered  in 
this  chapter.  A comparison  of  the  two  energy  spectrums  is  included  in 
the  discussion  of  the  results  in  Chapter  V.  It  should  be  noted  that 
throughout  this  discussion  the  term  kinetic  energy  implies  kinetic 
energy  per  unit  mass,  since  the  total  kinetic  energy  has  been  defined 
as 


KE 


2 


u 


(4.1) 


Equation  (3.6)  defined  the  time-averaged  kinetic  energy  of  the 
standing  eddies  (SKE)  for  each  grid  point  in  terms  of  i7  and  7.  How- 
ever, u and  V are  defined  by: 


u 


Uo  + u' 


(4.2) 


I 


and 


V = V + V' 
0 
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(4.3) 


where  u^,  are  the  time-and-longitudinally  averaged  quantities  defin- 
ed by  Equation  (3.5)  and  u',  v'  represent  deviations  from  the  average 
values.  Subsitutlon  for  u and  v In  Equation  (3.5)  yields 


SKE 


(uo  + u')^  + (vo  + v')^ 


(4.4) 


To  find  the  longitude-averaged  standing  kinetic  energy  for  a 
given  latitude.  Equation  (4.4)  needs  to  be  integrated  around  the  lati- 
tude circle.  Since  Uq  Is  Independent  of  longitude  and  u'  averaged 
around  a latitude  circle  Is  zero,  the  following  result  is  obtained  for 
the  u-component: 


- 2ir  —2  I ^.^2 

1 / * t!  (¥)■“ 


2ir 


(4.5) 


where  ^ represents  the  kinetic  energy  of  the  stationary  zonal  motion, 
wave  number  zero,  and  the  Integral  on  the  right  hand  side  gives  the 
total  longitude-averaged  kinetic  energy  of  the  standing  eddies. 

It  Is  possible  to  express  ir(xr)  In  terms  of  a Fourier  Series  as 
follows: 


u(xr)  = Aq  + y (A  cos  sin 


m=l 


— , 2irX|r'. 

SaP*  Tin--’ 


(4.6) 
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N = number  of  points  around  latitude  circle  (72  for  this 
analysis), 
n = N/2, 

and  r = index  denoting  the  point  on  the  latitude  circle. 

The  coefficients  and  8^,  for  m ^ o,  are  defined  as: 


T 2 - 2Trmr 

^ = N i N 

r=-n 


(4.7) 


and 


_ 2""^- 

Bm  " I 

r=-n 


2Trmr 


(4.8) 


Since  n is  even.  Bp  is  zero.  For  m = 0, 


n-1 


Aq  = J-  I Ur  and  Bp  = 0. 


r=-n 


Using  the  Fourier  Series  expansion  for  u and  Parseval's  Theorem, 
it  is  possible  to  express  the  longitude-averaged  standing  kinetic 
engergy  as: 


1 

k I (-r)"'  ■ 


1 -2 

^ “r 
r=-n 


2 ^m^l 


(4.9) 
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Comparing  Equation  (4.5)  and  (4.9)  leads  to  the  conclusion  that 


= u. 


(4.10) 


(4.11) 


From  Equation  (4.11)  It  Is  evident  that  the  contribution  to  the  longi- 
tude-averaged kinetic  energy  of  the  standing  eddies  for  a given  wave 
number  (SKE^,)  can  be  computed  as  follows: 


= 


(4.12) 


Longitudinal  averages  of  kinetic  energy  for  the  meridional  wind  com- 
ponents are  found  In  the  same  manner  by  replacing  u quantities  by  the 
corresponding  v quantities. 

Equation  (4.12)  was  used  to  calculate  the  values  of  SKE^,  for  wave 
numbers  0 to  36  at  fourteen  latitudes.  Values  for  u and  v components 
were  computed  separately.  The  results  for  wave  numbers  0 to  18  are  con- 
tained In  Table  3. 

Analysis  of  these  results  did  not  produce  anything  unexpected.  In 
mid-latitudes  much  of  the  standing  kinetic  energy  Is  contained  in  wave 
number  zero  of  the  u component.  At  200  mb  for  the  u component,  there 
is  a shift  of  the  energy  to  wave  numbers  2,  3,  and  4 near  the  equator. 

In  fact,  at  4°  S and  4°  N wave  number  zero  contains  essentially  no 
energy  at  200  mb.  However,  at  500  mb  the  maximum  energy  is  still  con- 
tained In  wave  number  zero.  The  Increase  In  the  total  200  mb  standing 
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kinetic  energy  at  4°  S is  a reflection  of  the  mid-Pacific  jet. 

For  all  latitudes  there  is  little  detectable  energy  above  wave 
number  5.  Another  contrast  between  200  mb  and  500  mb  is  the  latitude 
at  which  the  maximum  total  KE  occurs  in  the  Northern  Hemisphere.  At 
200  mb  the  maximum  occurs  at  28®  N latitude,  while  the  500  mb  maximum 
occurs  at  44°  N.  The  same  trend  appears  in  the  distribution  of  maxi- 
mum energy  for  wave  number  zero.  This  indicates  that  the  maximum 
kinetic  energy  at  200  mb  is  associated  with  the  sub-tropical  jet,  while 
the  maximum  energy  at  500  mb  occurs  in  conjunction  with  the  polar  jet. 
It  is  obvious  from  this  example  that  an  apparent  trend  at  one  level 
of  the  atmosphere  cannot  always  be  extrapolated  to  a different  level. 

The  V component  shows  a preference  for  wave  number  three  in  the 
Northern  Hemisphere  and  wave  number  four  in  the  Southern  Hemisphere  at 
both  200  mb  and  500  mb.  The  200  mb  total  column  shows  a maximum  in  the 
sub-tropics  in  both  hemispheres,  but  the  500  mb  v component  has  only 
one  significant  maximum  at  52°  N.  This  again  suggests  that  the  maximum 
energy  at  500  mb  is  associated  with  the  polar  front.  At  200  mb  it 
appears  there  is  about  as  much  standing  kinetic  energy  of  the  v compon- 
ent associated  with  the  polar  jet  as  with  the  sub-tropical  jet.  The 
results  obtained  from  the  spectral  analysis  of  the  DST  data  agree  with 
the  general  features  of  those  reported  in  other  investigations  (see 
i.e.,  Kao  and  Wendell,  1970). 

Also,  Winn-Nielson  (1967)  found  that  the  N.H.  maximum  of  the  zonal 
kinetic  energy  occurred  in  February  and  the  minimum  occurred  in  July. 
Thus,  it  can  be  assumed  that  the  values  given  in  Table  3 represents 
approximate  maximum  values  for  the  Northern  Hemisphere  and  minimum 
values  for  the  Southern  Hemisphere. 
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CHAPTER  V 

KINETIC  ENERGY  OF  TRANSIENT  WAVES 

The  time-and-longitude  averaged  kinetic  energy  for  the  transient 
waves  was  decomposed  in  wave  number  space  in  a manner  similiar  to  that 
described  in  Chapter  IV,  The  transient  part  of  the  flow  is  defined  as 
the  instantaneous  departure  from  the  time-averaged  flow  as  follows: 

u*  = u - IT  (5.1 ) 

and 

V*  = V - V (5.2) 

where  iT,  v are  defined  by  equations  (3.1)  and  (3.2), 

u,  V are  the  actual  wind  components  for  each  observation 
and  u*,  V*  are  the  components  of  the  transient  flow. 


Substitution  of  uj;  for  Uf«  in  Equation  (4.9),  results  in  the  follow- 
ing equation  for  the  longitude-averaged  transient  energy: 


B*  = I-  y u*  sin  Iimr 
^r=-n  N 


(5.5) 


For  m = 0,  A*  * i I u*  and  B*  = o 
” r=-  n 


An  quantities  except  the  starred  ones  are  defined  the  same  as  in 
Chapter  IV.  To  calculate  the  time-and-longitude  averaged  transient 
energy,  it  is  necessary  to  integrate  Equation  (5.3)  over  the  time 
period  in  the  following  manner: 


•I  T ^TT/  *v2  , I ll-l 

^ / / 2 ‘l^dt  = ^ / ( J -u*  )dt 


T n-1 


0 r=-n 


T / ( ^ * T I (AS^  * . 

0 m=l 


(5.6) 


Moving  the  integral  on  the  right  hand  side  of  Equation  (5.6)  inside 
the  wave  number  sunmiation  and  replacing  it  with  a summation  over  time 
yields: 


T , n-1 


M .*2 


n M 


(5.7) 


where  M = 70. 


If  it  is  considered  that  the  transient  flow  is  composed  of  a 


zonally  average  value  and  a deviation  from  this  average,  u*  can  be 
written  as 


u*  = u*  + u*' 

0 

where  u*  = transient  part  of  zonal  average 
and  u*'  = transient  eddy. 

Applying  the  above  definition  to  the  terms  on  the  right  hand  side  of 
Equation  (5.7)  results  in  A*^  = uj^  and 


1 

2irT 


T Zir,  -2 

/ / 

0 0 ^ 


'm=l  i=l 


Thus,  the  contribution  to  the  time-and-longitude  average  transient 
kinetic  energy  for  each  wave  number  (TKE„,)  can  be  computed  from 

ffij,  ■ <5.8) 

The  results  of  using  Equation  (5.8)  to  compute  the  values  of  TKE^ 
are  presented  in  Table  4.  TKe  latitudes  and  wave  numbers  are  the 
same  as  for  the  standing  eddies. 

Again  the  results  are  not  unexpected.  Table  4 shows  several 
contrasts  when  compared  with  Table  3.  The  total  transient  kinetic 
energy  of  the  u component  is  much  smaller  than  the  standing  kinetic 
enr-gy  in  mid-latitudes.  At  the  equator,  the  contributions  from  the  u 
components  are  about  the  same  magnitude  for  the  transient  and  standing 
kinetic  energies.  For  the  v components  the  transient  kinetic  energy  is 
significantly  larger  than  the  standing  kinetic  energy  at  all  latitudes. 


TABLE  4 

Decomposed  Time-averaged  Transient  Kinetic 
Energy  for  Wave  Numbers  0-18 
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Also,  the  distribution  of  energy  with  respect  to  wave  number  is  differ- 
ent for  the  transient  energy. 

The  transient  kinetic  energy  of  the  u component  is  concentrated 
in  wave  numbers  1 to  4 with  only  a small  amount  of  the  total  in  wave 
number  one.  For  the  v component  the  transient  kinetic  energy  is 
clustered  between  wave  numbers  3 and  8.  However,  in  the  tropics,  the 
transient  energy  is  spread  out  fairly  evenly  over  all  the  wave  numbers. 
This  implies  that  there  is  no  dominate  mechanism  in  the  tropics  in 
this  data  set  which  can  be  associated  with  the  transient  kinetic 
energy. 

It  should  also  be  noted  that  some  of  the  kinetic  energy  in  the 
larger  wave  numbers  may  be  due  to  gravity  type  waves.  Some  are  pro- 
bably real  but  others  may  be  induced  by  the  numerical  approximations 
of  the  model. 

The  supposition  that  the  model  is  producing  some  spurious  waves 
results  from  the  fact  that  at  some  latitudes  there  was  an  increase  in 
the  kinetic  energy  around  wave  numbers  33-36. 
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CHAPTER  VI 

SUMMARY  AND  CONCLUSIONS 

The  analyses  presented  here  Indicate  that  the  DST  data  assimilat- 
ed by  the  GISS  model  is  useful  for  global  studies  of  the  general 
atmospheric  circulation.  The  time-averaged  charts  presented  in  Chap- 
ter III  correlated  very  well  with  previously  produced  charts.  A close 
look  at  the  time-averaged  charts  reveals  that  this  data  set  is  only 
valid  between  about  60“  S and  60®  N.  Poleward  of  these  latitudes 
there  are  many  places  where  two  or  more  grid  points  have  the  same 
values.  Lack  of  sufficient  data  and  the  closeness  of  grid  points  pro- 
bably contributed  to  this  duplication.  Given  sufficient  data  to 
totally  initialize  the  grid,  the  GISS  model  should  provide  fine  qual- 
ity data  for  the  entire  globe.  A positive  feature  of  the  DST  data  is 
the  ease  with  which  it  can  be  handled  in  horizontal  arrays. 

Interpolc tion  of  the  wind  components  from  sigma  levels  to  constant 
pressure  levels  by  a cubic  spline  appears  to  have  been  very  successful. 
There  was  some  concern  as  to  the  effectiveness  of  interpolating  the  tt 
fields  to  wind  data  points,  however,  the  results  support  the  technique 
used  here. 

The  success  of  this  technique  suggests  it  could  be  used  to  gener- 
ate winds  at  most  any  level,  between  the  top  and  bottom  sigma  levels, 
for  which  winds  are  desired.  By  making  full  use  of  the  vertical  reso- 
lution of  the  model,  it  should  be  possible  to  see  the  time-and-longitude 
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averaged  meridional  circulation  in  more  detail. 

It  was  encouraging  to  see  significate,  though  small,  values  of 
the  time-and-longitude  averaged  meridional  wind  component  which  were 
organized  into  patterns.  Also,  the  time-averaged  values  of  the 
meridional  wind  show  organized  wave  patterns  in  both  hemispheres 
which  are  very  well  defined.  Thus,  the  winds  in  this  data  set  are 
much  less  constrained  than  in  earlier  studies. 

An  original  objective  of  this  research  was  to  study  the  DST 
divergence  fields.  The  time-averaged  divergence  showed  areas  of  organ- 
ization, especially  at  200  mb,  indicating  that  the  DST  data  can  be 
used  for  divergence  studies.  However,  it  is  uncertain  how  much  detail 
can  be  resolved  in  the  divergence  fields.  The  randomness  observed  in 
some  areas  of  Figures  26  to  31  may  be  due  to  some  partial  contamin- 
ation by  gravity  type  waves.  This  conclusion  results  from  the  fact 
that  an  increase  was  observed  in  the  kinetic  energy  for  high  wave 
numbers  at  some  latitudes.  Another  source  of  uncertainty  for  both 
the  vorticity  and  divergence  fields  is  the  mountain  effect  on  the 
fields.  This  effect  was  discussed  in  Chapter  II  and  was  probably 
greatly  reduced  by  the  interpolation  to  constant  pressure  levels. 

Wave  number  decomposition  of  the  kinetic  energy  spectrum  around 


latitude  circles  is  easy  to  accomplish  because  of  the  structure  of  the 
grid.  Although  there  were  no  surprises  in  the  results  from  Chapter  IV 
and  Chapter  V,  it  was  noted  that  there  is  some  difference  between  the 
relative  amounts  of  energy  carried  by  waves  of  a certain  wave  number 
at  different  levels.  Thus,  to  adequately  describe  the  energy  trans- 
port in  the  atmosphere,  it  is  important  to  consider  several  vertical 
levels. 
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Most  of  the  features  observed  in  the  data  conformed  to  previously 
calculated  averages.  The  200  mb  mid-Pacific  trough  over  the  equator 
stands  out  as  a subject  worthy  of  future  research.  From  this  data  it 
appears  that  flow  from  south  of  the  equator  is  feeding  into  the  sub- 
tropical jet  which  extends  across  Central  America.  A study  of  the 
energy  transfer  involved  in  this  process  could  prove  to  be  very  inter- 
esting and  informative. 
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